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Abstract. By diagonalization of a generalized superoperator for solving the master equation, we investi- 
gated effects of dissipative and dephasing environments on quantum state transfer, as well as entanglement 
distribution and creation in spin networks. Our results revealed that under the condition of the same de- 
coherence rate 7, the detrimental effects of the dissipative environment are more severe than that of the 
dephasing environment. Beside this, the critical time tc at which the transfer fidelity and the concurrence 
attain their maxima arrives at the asymptotic value to = n/2X quickly as the spin chain length A'^ increases. 
The transfer fidelity of an excitation at time to is independent of A'' when the system subjects to dissipative 
environment, while it decreases as A'^ increases when the system subjects to dephasing environment. The 
average fidelity displays three different patterns corresponding to A'^ = 4r -I- 1, A'^ = 4r — 1 and TV — 2r. 
For each pattern, the average fidelity at time to is independent of r when the system subjects to dissipa- 
tive environment, and decreases as r increases when the system subjects to dephasing environment. The 
maximum concurrence also decreases as N increases, and when A' — >■ 00, it arrives at an asymptotic value 
determined by the decoherence rate 7 and the structure of the spin network. 

PACS. 03.67.-a Quantum information - 03.67.Mn Entanglement production, characterization, and ma- 
nipulation - 03.65.Yz Decoherence; open systems; quantum statistical methods 
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1 Introduction 

' An important and new emerging pursuit of quantum infor- 
] mation processing (QIP) task is the high-fidelity transmis- 

■ sion of quantum states across a pre-engineered quantum 
' spin network. This is so because when performing any 

QIP tasks, one needs to exchange quantum information 

■ between distant nodes (e.g., the core processor, storage, 
\ etc.) of a quantum computer. In general, the basic idea 

■ of transmitting a quantum state from one location of a 
' network to another proceeds in three steps [1]. The first 
, step is the initialization of the spin medium in a fiducial 

■ pure state, then the sender Alice encodes the state needs 
\ to be transmitted at one node of the network, and time 

evolution of the system for a proper interval of time, fi- 
nally, the state will be recovered at the nodes belonging 
to the receiver Bob with certain fidelity. Particularly, the 
quantum channel may enable state transfer with fidelity 
better than any classical communication channel 

Since the seminal work of Bose 1 , in which the author 
demonstrated that an unmodulated ferromagnetic spin 
chain with Heisenberg interactions can be used as a chan- 
nel for short distance quantum communication, there has 
been an increasing discussion of quantum state transfer 
in various physical systems in the literature. Among such 
systems, the solid-state quantum system with spin-spin in- 
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teractions offers an excellent theoretical framework for de- 
signing networks for perfect state transfer (PST) [1,2,31H1 
RR[71R[^[TT)l[TTl[ni[ni[ni[TK|Pi)] ■ Particularly, Christandl 
et al. suggested a PST algorithm which can transfer an ar- 
bitrary quantum state between the opposite ends of a spin 
chain or the two antipodes of the one-link and the two- 
link hypercubes Zhang and Long et al. [13] realized 
this PST algorithm in a three-qubit XX spin chain using 
liquid nuclear magnetic resonance (NMR) system. Since 
then, quantum spin networks as an ideal communication 
channel to realize QIP tasks have been extensively stud- 
ied. Shi et al. |6| presented a series of perfect spin chan- 
nels according to the spectrum-parity-matching condition 
they derived. Then, Kostak et al. [2] established a gen- 
eral formalism for engineering spin Hamiltonians for PST 
in spin networks of arbitrary topology and coupling config- 
uration. Christandl's innovative works OH] were extended 
by Jafarizadeh and Sufiani in a recent work |15| , in which 
they adopted distance-regular graphs as spin networks and 
found that any such network can achieve unit fidelity of 
state transfer over arbitrarily long distances. Moreover, 
D'Amico et al. [TU] showed that one can create and dis- 
tribute entanglement with an interaction-modulated Y- 
shaped spin network, particularly, with a slightly compli- 
cated bifurcation structure, one can even froze the entan- 
glement by applying a phase flip to one spin out of each 
pair. 

Apart from the aforementioned protocols which mainly 
concentrated on quantum spin chains with only nearest- 
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neighbor (NN) couplings, in Ref. [17^ Paternostro et al. 
presented another scheme which enables nearly optimal 
state transfer in imperfect artificial spin networks with all 
the qubits are mutually coupled. Then Kay demonstrated 
that PST is also possible in the presence of next-nearest- 
neighbor (NNN) couplings [T5]. Moreover, compared to 
the case where the system contains only two-spin inter- 
actions, the authors in Ref. [T^ showed that the speed 
of quantum state transfer across the XY spin chain can 
be significantly increased by introducing the three-spin in- 
teraction. Besides the spin-1/2 systems, state transfer ef- 
ficiency of a bilinear-biquadratic (BB) spin-1 Heisenberg 
chain has also been discussed [20] . 

Most recently. Franco et al. [21] presented a control- 
limited scheme [35] for PST through a pre-engineered spin 
chain. Different from the previous schemes whose achieve- 
ments relies crucially on the initialization of the channel 
state. Franco et al. [H] demonstrated that such state ini- 
tialization is in fact inessential to the performance of the 
protocol if proper encoding at the end of the chain is per- 
formed. The key requirements for their scheme are the ar- 
rangement of proper time evolution and the performance 
of clean projective measurements on the two end spins, 
which considerably relaxes the prerequisites for obtaining 
reliable state transfer across interacting-spin systems. Mo- 
tivated by this innovative work, Markiewicz and Wiesniak 
[22] proposed a special type of encoding strategy for PST 
of a two-qubit state, where no remote-cooperated global 
state initialization and any additional communication are 
needed. 

Besides these exciting progresses in this direction, how- 
ever, one cannot neglect the fact that real quantum sys- 
tem is very susceptible to decoherence 24,,25 , 26. 27, 28, 29] 
l30U31ll32ll33ll34ll35j . For nearly all QIP tasks, the quan- 
tum system of interest will inevitably be interacted with 
its surrounding environments (or the thermal reservoirs). 
This unavoidable mutual interaction often cause the initial 
state of the system becomes entangled with the environ- 
ments in an uncontrollable way, and it is this entanglement 
of the system with the environment that induces decoher- 
ence. The decoherence can affects quantum interferences 
of quantum systems and leads to the degradation of quan- 
tum coherence, and thus becomes a serious limiting factor 
bafffing the physical realization of a quantum computer. 
Since the system of interest is inevitably subject to de- 
coherence and decay processes, no matter how much they 
may be screened from the external environments, it would 
seem important to consider possible methods to minimize, 
delay, or even eliminate this unwanted detrimental effects 
of environments in the practical realization of various QIP 
tasks I3B1I37]. 

Dissipation and dephasing are two important sources 
of decoherence in practice, and their effects on entan- 
glement dynamics have been analyzed very recently in 
Refs. [56,27,28,22] for the special case of two spin-halves 
where an analytical treatment is possible. The results re- 
vealed that the entanglement dynamics depends not only 
on the parameters of the system, but also on the system- 
environment coupling strength and the initial states of 



concern. How these decoherence modes affect quantum 
state transfer fidelity in a spin chain is an interesting and 
urgent problem needs to be solved, and it is also of both 
theoretical and experimental significance. We will address 
this question in the present paper. The contents are or- 
ganized as follows. In Section 2, we introduced the tech- 
nique for solving the Lindblad form of the master equa- 
tion. Based on this technique, we obtained the time evo- 
lution operator of the system and the explicit form of the 
state transfer fidelity. Then in Section 3, we investigated 
dynamics of the state transfer fidelity with the system 
subjects to dissipative and dephasing environments, where 
we concentrated on several ideal spin channels in the ab- 
sence of decoherence, and aimed at revealing the extent 
to which the decoherence affects quantum state transfer. 
Section 4 is devoted to the analysis of the effects of the 
dissipative and dephasing environments on entanglement 
distribution and creation between distant parties through 
a pre-engineered spin network. Finally, we concluded this 
paper with a short discussion in Section 5. 



2 Solution of the master equation 

As pointed in Section 1, the unavoidable mutual interac- 
tion of an open quantum system with its surrounding en- 
vironment is an important source of decoherence. In order 
to describe such phenomena, a master equation approach 
can be used. In the present work, we investigate how the 
dissipative and the dephasing environments affect quan- 
tum state transfer, as well as entanglement distribution 
and creation through a pre-engineered spin network. We 
assume that during the decoherence processes each spin of 
the system interacts only, and independently, with its own 
environment (this assumption is legitimate provided the 
constituents composing the quantum system are separated 
by distances large enough, however, it may not be totally 
valid in the case where the spins are close enough to each 
other to couple strongly enough), then in the Markovian 
limit, the dynamics of the system can be described by a 
general quantum master equation (ME) of the Lindblad 
form [SHllSni 

f^^-^[H,p] + J2Cr.p, (1) 

n 

where H and p denote, respectively, the Hamiltonian and 
the density operator of the system, and the summation 
on the second term of the right-hand side of this equation 
runs over all the spins involved. refers to the Lindblad 
superoperator, which describes the independent interac- 
tion between each spin and the environment and can be 
written, under weak system-environment interaction and 
in the Markovian limit, as ^B] 

^nP = ^(2c„pct - clcnP - pclcn ) . (2) 

Here 7 is the phenomenological parameter that describes 
the coupling strength of a qubit with its local environment. 
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The system-environment coupling operator is given by 
Cn = (Jn for the dissipative environment, and c„ — cr^cr" 
for the dephasing environment. The first term on the right- 
hand side of Eq. (1) generates a coherent unitary evolu- 
tion of the system, while the second term represents the 
decoherence effects of the environment on the system and 
generates an incoherent dynamics of the system. 

In order to assess the extent to which the decoherence 
affects dynamics of a correlated spin system, one must 
first try to obtain solutions of the ME. Although for sys- 
tems with certain symmetries, the ME can be solved by 
adopting the superoperator method [ID] or the Lie alge- 
braic method (see, e.g., Ref. [H] and references therein). 
For general cases, however, this is still an uneasy task 
because the ME expressed in Eq. (1) contains terms cou- 
pled with each other in a non-trivial way. Although one 
can convert it to an equivalent c-number partial differen- 
tial equation [42^ and then solve it numerically using the 
conventional Runge-kutta algorithm; however, for system 
with large number of spins, this is still a troublesome task 
for it wasting too much computation resources to obtain 
the density operator p{t), and the situation becomes even 
more cumbersome when one performs simulations of p{t) 
with different initial states because for every initial state 
one must perform the Runge-kutta algorithm again to ob- 
tain the solutions of the ME. In order to circumvent this 
quandary, in the present paper we resort to another tech- 
nique to solve the ME, as one will see, the success of this 
technique depends solely on the exact diagonalization of 
a generalized superoperator A, and it has no relation with 
the initial state of the system. This technique provides an 
efficient way for the investigation of the dynamics of the 
system and not only allow for an efficient numerical anal- 
ysis, but also, in some cases, for exact algebraic solutions. 

The basic idea of this technique is as follows. Instead of 
solving the ME directly, one can first define two ancillary 
operators A and p, which satisfy the following relations 



can be divided into the following three terms, i.e. 



where 



dp 



P = (Pl,:,P2,:, • • ■■.Pd-.f 



(3) 



(4) 



Here d denotes the dimension of the density matrix p, 
(. . -Y" denotes the transpose of (. . .), and pn.- — [pni, Pn2, 
. . . , pnd\ denotes the elements lie in the nth row of p. After 
defining these two ancillary operators, the formal solution 
of Eq. (3) can be expressed explicitly as 



p{t) = u{t)m- 



(5) 



Since p is a d x d density matrix, the time evolution op- 
erator U{t) = exp(ylt) is a d^ X matrix, and the initial 
state of the system p(0) is a d^ x 1 column vector. More- 
over, from the definition of p{t) one can obtain the equality 
P{i-i)d+jit) = Pij{t). 

Now the key problem is how to obtain the explicit 
expression of A. The crucial observation to this problem 
is that the right-hand side of the ME expressed in Eq. (1) 



(6) 



where A and B are operators which can be obtained di- 
rectly from Eqs. (1) and (2) as A — —iH — ^^^{cl^Cn/2), 

B = iH — 7X]n(CnCn/2). Hereafter we adopt the conven- 
tional notation Mij to denote the element lies in the ith 
row and jth column of a matrix M, and Vi denotes the ith 
element of a column vector V. In order to express the ME 
in forms of dp/dt = Ap, we make the following three trans- 
formations: Ap A^^^p, pB ^ A^'^^p and 7X)nCnPcJj >-)■ 

A^'^^p. To ensure the achievement of the first transforma- 
tion Ap i-> A^^'^p, the equality {Ap)ij — 

must be hold, or equivalently, J2k Ai,kPk,j = 

then after a straightforward algebra, one can obtain yl^^^ = 
A® I . Similarly, for the second transformation, one can 
obtain yl^^^ = I ® B' , while for the third transformation, 
one can obtain A^^^ ~ ® ^n- the above expres- 

sions, B' denotes the transposed matrix of B, c* denotes 
the conjugate of c„, and / is the d x d identity matrix. 
With the help of these results, Eq. (3) becomes 



dp 



(7) 



This is just the desired form of the transformed master 
equation, where the newly defined operator is given by 
A = A®I + I®B' + ^Y.n ® <• Once the exphcit 
form of A has been found, the time evolution operator 
U{t) = exp {At) can be formed and the formal solution of 
the master equation is obtained, which can be expressed as 
Eq. (5). To obtain the density operator at arbitrary time t, 
one only needs to apply the time evolution operator U (t) 
to the initial state p(0) of the system. 

Now we try to convert the matrix exponential form of 
the time evolution operator U (t) — exp (At) to the matrix 
form. For this purpose, we expanded it in terms of the 
Taylor series as 



exp (At) = X 



(Aty 



(8) 



where n\ denotes the factorial of n. If the d^ x d^ matrix A 
is diagonalizable {A cannot always be diagonalized since 
it is generally complex and non-Hermitian), i.e., there ex- 
ists an invertible matrix M that satisfying the relation 
M~^AM = diag{£i,£2, ■ ■ ■ , £d'^)j with the diagonal entries 
£k (fc = I7 2, . . . , d^) being the eigenvalues of the operator 
A, then from Eq. (8) one can obtain 



exp {At)M = 



(M-^AtMY 



= cxp{£t). (9) 



From Eq. (9) one can obtain directly that the matrix form 
of the time evolution operator U (t) is given by 



ij{t) = Mexp{£t)M-^, 



(10) 
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where represents the inverse matrix of M, and M is 
the matrix with the eigenvectors of A as its columns. 

From the above analysis, one can see that the achieve- 
ments of this technique depends solely on the exact diag- 
onalization of the operator A, once A is diagonalized, one 
can obtain p{t) as well as the matrix form of p{t) at any 
time t by applying the time evolution operator U{t), par- 
ticularly, U{t) has no relation with the initial state p(0) 
of the system, which is economical for performing numer- 
ical calculations with different initial states, and this is 
different from that of the Runge-kutta algorithm. 

In the present paper, we first consider effects of dissi- 
pative and dephasing environments on state transfer in a 
spin chain. We assume the quantum state to be transmit- 
ted is encoded at the mth spin as \(pin) — cos (0/2) |0) -I- 
e*"^ sin (0/2)|l) (with |0) and |1) represent the state of spin 
up and down, respectively), and all the other spins in the 
chain are initialized to the ground state |0), thus the initial 
state of the whole system at time t = becomes 



IV'(O)) =cos-|0) 



e**^ sin - I to) 



(11) 



where |0) = |00...0), |to) = cr+ 10)®^, ^ 6i tt and 
^ (p ^ 27T are the relative phase angles. For this type of 
initial state, the nonzero elements of p(0) can be written 
as 

PNm~N+27n-l{0) = siu^ f , PN-^+2N+liO) COS^ f , 

(12) 

Except these four elements, all the other elements of p{0) 
are zero. Here A'' denotes the length of the chain, and 
we have used the fact that for initial state \ijj{0)) with 
the system subjects to dissipative and dephasing environ- 
ments, its dynamics is completely determined by the time 
evolution in the zero and single excitation subspace "Hoei 
PBl.27.28,29 . thus it suffices to restrict our attention to 
the dynamics of p{t) in this {N + l)-dimensional subspace 
spanned by {|0), |1), |2), . . . , \N)} which greatly facilitates 
the following computation process. 

To evaluate the extent to which the decoherence af- 
fects state transfer in a spin chain, we adopt the concept 
of fidelity 1 / = {yyin\Pnit)\'Pin) as an estimation of the 
quality of the state transfer from the sender, convention- 
ally named Alice, to the receiver Bob, where Pn{t) is the 
single qubit reduced density matrix. The fidelity measures 
the overlap between the input state and the output state. 
In general, if one encodes the state needs to be transmitted 
at one end of the chain, then after some time t, the state 
will be recovered at another end of the chain automati- 
cally, with however, the fidelity / < 1. With elaborately 
designed structures of the spin chain which is ideally pro- 
tected from its surrounding environments, the maximum 
fidelity may reach unity pHllSl fTlilTS] . When subjects to 
decoherence environments, as can be seen in the following 
sections, the transfer fidelity cannot reach unity even with 
elaborately designed spin structures. 



In the standard basis {|0),|1)}, the single qubit re- 
duced density matrix Pn{t) can be obtained by tracing 
out all other qubits except n from p(t) as 



/.N 1 - Pnn{t) PN+I,n{t) 

P''^^>^ \pn,N+l{t) Pnn{t) 



(13) 



Here Pmn{t) denotes the element lies in the TOth row and 
nth column of p{t). Then the state transfer fidelity from 
the TOth qubit to the nth qubit can be computed as 



/ = COS^ - - PnnC0s9+ -{e"^PN+l,n + 6 Pn.N +l) Sin . 

(14) 

In terms of p{t) the elements of Pn{t) can be written as 

Pnn{t) = PNn-N+2n-l{i)-, 

Pn,N+l{t) = PNn+n{t), 
PN+l,n{t) ^ PN^+N+n{t)- (15) 



The exphcit form of pNn-N+2n-i{t), PNn+n{t) and 
PN^+N+n{t) can be expressed as 



PNn-N+2n-l — Unu- N+2n-l,Nm- N+2m-l Sm — 

fr 2 9 

+ U Nn-N+2n-l,N'^+2N+l COS - 

+ -^{UNn~N+2n-lMm+me^'^ 

+ UNn-N+2n-l,N^+N+m^ *"^)sin0, 
■2^ 

PNn+n — UNn+n,Nm,-N+2m-lSill — 

2 ^ 

+ U Nn+n,N^+2N+l COS - 

e '*')sin6', 

■ 2^ 

PJV2+iV+n — t^Ar2+Af+n,Arm-A'+2m-l sin - 

2 ^ 

+ UN^+N+n,N^+2N+lC0S - 



+ UN^+N+n,N^+N+me *'^) siu 



(16) 



Substituting these results into Eq. (14), one can obtain 
the explicit form of the state transfer fidelity. Particularly, 
when 9 = TT, we obtain the transfer fidelity of an excita- 
tion as / = (l|p„(t)|l) = UNn-N+2n~l,Nm~N+2m-l- On 

the other hand, in order to assess the efficiency of the 
quantum spin channel of interest, it is more beneficial to 
calculate the average fidelity (the fidelity averaged over all 
pure input states in the Bloch sphere) F = J fdfl = 
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lo lo^ OdOdcf) [T] , which can be computed as 



3 State transfer in decoherence spin channels 



F 



UN2+N+n,N^+N+m + UNn+n,Nm+m 



+ 



6 



6 

UNn-N+2n-l,N^+2N+l 1 

6 2' 



(17) 



From the above equation one can see that the average 
fidehty F is solely determined by the time evolution op- 
erator U (t) . In a previous work [33] , Bowdrey et al. have 
also derived a similar expression for average fidelity in 
the form of unitary (or anti-unitary) operators. Moreover, 
Eq. (17) can be simplified under certain special condi- 
tions. For example, if one chooses the correlated spin chain 
Hamiltonian of the considered quantum system as H = 

then it is direct to show that the last column of H can be 
written as H- d = (0, 0, . . . , 0, c)-'" in the standard basis, 

where d is the dimension oi H, c = X]ri(^n,"-i-i + ^n) 
is a nonzero number. Combination of these results with 
the expression of the operator A — A®I + I®B' + 
7 X]„ c„ (g) c* , one can show that the last column of A 
has the form A. j2 = (0, 0, ... , 0)-^, thus by adopting the 
Taylor series expansion of the time evolution operator 
U{t) = exp(yli), one can compute the last column of U{t) 
as U-^cP = (0; 0, . . . , 0, 1)^ (note that this expression is 
valid for all subspaces, when restricted to the subspace 
7^001 J it holds true even for the system Hamiltonian with 
inhomogeneous x and y components), this gives rise to a 
simplification of Eq. (17) as 



F = 



6 

, UNn-N+2n-l,Nm-N+2m-l 1 

+ Q +2- (18) 



Eqs. (17) and (18) are the main results we obtained 
in this paper, which will be used in the latter discus- 
sion of quantum state transfer when an explicit form of 
the spin chain Hamiltonian is given. Particularly, for the 
mirror-symmetric Hamiltonian Hm (by mirror symmetry, 
we mean that the interaction Hamiltonian of the spin 
chain has symmetric coupling strengths about the cen- 
tre qubit or the centre link, i.e., the coupling strength 

N-n+l) we have C/jV^+Ar+n^A'^ + TV+m = 



In the absence of decoherence environments 



Nn+n^Nm+m ' 

(i.e., 7 = 0), this yields UM^+M+n,N'^+N+m + UNn+n,Nm+m 

= 2Re{(n| exp(— i77mi)|"^)} = 2|(n| exp(— zi?TOt)|m) | cos a 

and UNn-N+2n-l,Niyi-N+2m-l = | (tI | exp(-ii/„i) |to) | ^ , 

where a = arg{(?i| exp(— iiJ,„t)|m)} denotes the argument 
of the complex number {n\eyi^{—iHrnt)\m). Clearly, this 
is just that of Eq. (6) in Ref. [T], which describes average 
fidelity in the non-disturbed case. 



Recently, intense research efforts has been devoted to the 
interacting spin systems which were proposed as poten- 
tial candidates to simulate the relation between qubits in 
a quantum computer 03]. This choice is due to the fact 
that such systems can be easily manipulated (e.g., by tun- 
neling potentials or energy bias) and scaled up to large 
registers. In the present paper, we investigate state trans- 
fer in this kind of system, with the addition of the pres- 
ence of the dissipative and dephasing environments. We 
concentrated on the spin channel with fixed but different 
coupling strengths between neighboring qubits, with the 
Hamiltonian given by 



( „x „x 



(19) 



n=l 



where cr" (a = x, y, z) are the usual Pauli operators acting 
on the nth spin, J„,„+i = \^Jn{N — n) is the modulated 
exchange coupling strength, and A is a scaling constant. 
This Hamiltonian is identical to the representation of the 
Hamiltonian of a fictitious spin S ~ [N — l)/2 particle: 
Hs = XSx, where Sx is its angular momentum operator 
in the x-direction. 

For this modulated spin chain, it has been shown that 
one can achieve perfect state transfer between the input 
node n and the output node TV — n -I- 1 after a time 
^0 — 7r/2A and at intervals of tt/A thereafter in the absence 
of decoherence environment, i.e., when t = {2k — l)7r/2A, 
k = {1,2,...}, the transfer fidelity reaches unity [3113] ■ 
When the decoherence is present, however, this ideal com- 
munication channel may be destroyed. To show that this 
is so, we consider effects of dissipative and dephasing en- 
vironments on quantum state transfer from one end of the 
chain to another. The average fidelity can be obtained di- 
rectly from Eq. (18) by replacing (m, n) with (1, N). More- 
over, since the Hamiltonian H expressed in Eq. (19) is mir- 
ror symmetric, we have UN^+2N.m+N+i = ^ir^+^r^^v+i' 
thus the average fidelity can be rewritten as 

F= — + + -. (20) 

Physically, one can maximize the average fidelity F by 
performing a phase flip operation or applying an external 
magnetic field along the z axis such that a is a multiple of 
27r. Furthermore, for the case of dissipative environment 
one can show that Ujs[2+m-i,i = \Un^+2N,n^+n+i\'^ , thus 
the expression of the average fidelity F is very similar to 
that of Eq. (6) in Ref. [1], which describes the average 
fidelity in the non-disturbed case. However, for the case 
of dephasing environment, one cannot obtain such a sim- 
ple relation. Except this, there exists another difference 
between the dissipative and the dephasing environments 
which may be useful in the following discussion. To demon- 
strate this explicitly, it is illuminating to rewrite the an- 
cillary operator A = A®I + I®B'+^ ^„ c„ ® c* as 
A = A{1) + yl(2), where A{1) ^ -i{H ® I - I ® H') and 
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^(2) =7E„{cn®c* -i[(4c„)(»/ + /®(4c„)']}, then 
for the case of dissipative environment, one can show di- 
rectly that A{1) and ^(2) commutes with each other, i.e., 
yl(2)] = 0, while for the case of dephasing environ- 
ment, one still cannot obtain this commutation relation. 

As a heuristic analysis, we first consider dissipative ef- 
fects on fidelity of quantum state transfer from one end 
of the chain to another at time to = {2k — 1)tt/2X. This 
choice of time to is based on the fact that for weak system- 
environment coupling strength 7, the difference between 
to and the critical time tc at which the transfer fidelity at- 
tains its maximum value is minimal (e.g., for N = 2 and 
7 = 0.1, the deviation for transfer fidelity of an excitation 
at time to is about 0.2496% from that at time tc), partic- 
ularly, for large N and small 7, this difference can even be 
neglected (see the inset of Fig. 1). For this kind of decoher- 
ence mode, from the commutation relation [yl(l), yl(2)] — 
one can find that the time evolution operator U{t) = 
exp(ylt) can be rewritten as U{t) = exp[yl(l)t] exp[yl(2)t]. 
For the first term, from the equality {A (S) L){B (g) M) = 
AB (g) LM and the Taylor series expansion of it one can 
obtain exp[yl(l)t] = [exp(-«t7f) (g)/][J(g)exp(it7f')]. When 
restricted to the zero and single excitation subspace T-Lqqi, 
the eigenvalues and the corresponding eigenvectors of H 
can be readily obtained as pS] 



£0 - 0, £k 



-(iV-2fc + l)A, 



1-00) = |0), \il>k) = T.n=i Cfc,«|n-), 



(21) 



where k = {1, 2, . . . , N}, and the coefhcient Ck,n is given 
by the following recursion relations 



1 



Cl,l 



V2(^-i) ■ 



£kCk,n-i - J{n-2){N -n + 2)ck.n~2 , ^ 
Cfc,« = , , ■ \n ^ 2). 



V(n-l)(VV-n + l) 



(22) 



Then exp[yl(l)to] can be obtained analytically (for the 
ease of presentation, we omitted its explicit expression 
since this provides no new insights on the system), and 
the (N'^ + N - l)-th row of it is [1, 0, 0, . . . , 0], while the 
(iV2 + 2iV)-th row of it is [0, . . . , 0, i(^'-i)(2fe-i) , 0, . . . , 0] , 
with the only nonzero element locating at the 

(iV^ + N + l)-th column. 

On the other hand, the eigenvalues and the eigenvec- 
tors of M2) can also be obtained analytically, the results 
are expressed as follows 

£; = 7(0, ci,ci,...,ci,C2), 
/ 3 



M 



\l d d 



\ 



9 

dj 



(23) 



where the elements of the (Af + l)^-dimensional row vector 
E are eigenvalues of i(2), M is the [{N + if x (A^ + 



l)^]-dimensional matrix with eigenvectors of A{2) as its 
columns. Except the parameters 1, d and g, all the other 
elements in the matrix M are zero. The explicit form of 
ci , C2 , d and g are given by 

ci = -(1,1,..., 1,0.5), C2 = -( 0.5,0.5 ■,0.5 ), 

N N 

d = (-1/72, 0, 0, . . . , 0), g = diag{l/V2, 1, 1, . • . , 1). 



JV-l-l 



N+l 



(24) 



Thus the analytical form of exp[yl(2)to] can also be ob- 
tained. The first and the {N^ + N + l)-th column of it 
can be written explicitly as [e"'''*", 0, . . . , 0, 1 — e"'*'*"]"^ and 
[0, . . . , 0, e-T*«/2, 0, . . . , 0]^, with the only nonzero element 
Q-ito/2 locating at the {N^ + N + l)-th row. 

From the above analytical results, one can obtain the 
transfer fidelity of an excitation and the average fidelity 
at time to = {2k — l)7r/2A as 



/(to) = e-^*", 

e^T'o/^cosa 
^ (*o) - 5 



6 



(25) 



where cos a = Re[i'^^~^^*-^'^~"'^^]. From the above equation 
one can see that when the system subjects to dissipative 
environment, the transfer fidelity /(to) of an excitation is 
independent of the spin chain length N. For fixed to, it 
will exponentially decay with the increase of the decoher- 
ence rate 7, while for fixed 7, it will exponentially decay 
with the increase of to (this phenomenon may be under- 
stood physically from that fact that we have assumed the 
dissipation acts independently and equally on each spin, 
and since there is a single excitation, the survical probabil- 
ity of this excitation should be an exponential function of 
the time to). The average fidelity F(to) is, however, depen- 
dent on the spin chain length N due to the existence of the 
phase factor cos a. As can be found from Eq. (25), F{to) 
displays three different dynamical patterns corresponding 
to A^ = 4r + 1, A^ = 4r - 1 and Af = 2r (r = 1, 2, . . .). 
For the latter two patterns of A^, one needs to perform 
a phase flip operation at the end spin belonging to the 
receiver [4] or apply an external magnetic field to spins 
of the system [1] to maximize the average fidelity (after 
these operations, the average fidelity at time to is also 
independent of the spin chain length A^ and will display 
behaviors completely the same as that depicted by the red 
diamonds in the bottom panel of Fig. 1). If there are no 
additional operations are performed, from Fig. 1 one can 
observe that only when A^ = 4r -I- 1 (or equivalently, A^ — 1 
to be divisible by 4) can the quantum spin channel may be 
superior to its classical counterpart, for which the highest 
fidelity for transmission of a quantum state is 2/3 [IB] . 

For arbitrary evolution time t, after a tedious compu- 
tation, we obtain the state transfer fidelity of an excitation 
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Fig. 1. (Color online) Transfer fidelity of an excitation (a) 
and the average fidelity (b) versus the rescaled time to/n = 
(2fe— 1) /2 (fc = 1, 2, . . .) with the system subjects to dissipative 
environment, where the decoherence rate and the scaling con- 
stant are given by 7 = 0.1 and A = 1. In (b), the plots from top 
to bottom correspond to A'' = 4r + 1, A*" = 2r and = 4r — 1 
(r = 1,2,...), with the straight line at _F = 2/3 shows the 
highest fidelity for classical transmission of a quantum state. 
The inset in (a) shows tc/n versus TV, tc is the critical time at 
which the transfer fidelity of an excitation (blue circles) and 
the average fidelity (red plus signs) attain their maxima. 



and the average fidelity analytically as 

/(i) = e-^*sin2(^-i)At, 

'^-\Xt)cosa e-^*sin2(^^i)Ai 1 



Fit) = 



6 



(26) 



From the above equation one can see that under the 
detrimental influence of dissipative environment, both f{t) 
and F{t) behave as suppressed oscillations as the time 
evolves. Moreover, the decay of the average fidelity F{t) 
displays three different patterns mediated by a phase fac- 
tor cosa = Re[i(^-i)(2fc-i)], 

From Eq. (26) one can also see that there exists a re- 
gion in which the dissipative spin channel is superior to its 
classical counterpart. This region can be obtained analyti- 
cally by solving the nonlinear equation F{tc) > 2/3, which 
gives rise to 7 < ln(2 cos^ a + 2 cos a-y/l -I- cos^ a + l)/tc- 
If we choose cos a = 1 and tc — to = n/2 (i.e., A = 1), 
then we have 7 < 1.122, and this result is independent of 
the spin chain length N. 

When the system subjects to dephasing environment, 
however, one cannot obtain an analytical result as that of 
the dissipative environment because [yl(l), yl(2)] ^ 0, so 
we resort to numerical techniques. The typical results are 
shown in Fig. 2, from which one can see that the transfer 
fidelity of an excitation f{to) also decays with the increase 
of to, however, when to — ^ 00, /(to) arrives at a nonzero 
value f{to — >■ 00) = 1/N, and this is different from that of 
the dissipative environment, in which /(to) decays to zero 
asymptotically. There is another difference exists, namely, 
for any fixed decoherence rate 7, /(to) decreases with the 
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Fig. 2. (Color online) Transfer fidelity of an excitation (a) 
and the average fidelity (b) versus the rescaled time to/n = 
{2k — l)/2 (fc = 1, 2, . . .) with the system subjects to dephas- 
ing environment. Here the decoherence rate and the scaling 
constant are given by 7 = 0.1 and A = 1. For (a), the plots 
from top to bottom correspond to TV = 2 to TV = 20. For (b), 
the blue diamonds from top to bottom correspond to r = 1 
to r = 4, the black crosses from top to bottom correspond to 
r = 1 to r = 10, and the red squares from top to bottom cor- 
respond to r = 1 to r = 5. The inset in (b) is an exemplified 
figure plotted for the comparison of the average fidelity F be- 
fore (hollow squares) and after (solid squares) the phase shift 
caused by the factor cos a being corrected. 



increase of TV for the case of dephasing environment, which 
is disadvantageous for long distance quantum communica- 
tion. The average fidelity F(to) is also dependent on the 
spin chain length, and it displays three different patterns 
for different kinds of TV. As can be seen from Fig. 2b, only 
for the case of TV = 4r + 1 (r = 1, 2, . . .) can the quantum 
spin channel may be superior to its classical counterpart 
if no other performance have been made. When to — > 00, 
the average fidelity also arrives at a constant value which 
is also dependent on TV, and can be computed analytically 
as F(to ^oo) = l /6TV +1/2. 

For fixed spin chain length TV and critical time t^, there 
also exists a critical decoherence rate 7c before which the 
dephasing channel is superior to any classical communi- 
cation channel, for which can only ensure transmission of 
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Fig. 3. Critical decoherence rate 7c versus A'^ with the system 
subjects to dephasing environment, where the scahng constant 
A has been chosen to be 1, and the phase shift for N = Ar — 1 
and N — 2r has been assumed to be ehminated. 
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Fig. 4. (Color online) FWHM At of the first oscillation versus 
A'^ with 7 = 0.1 and A = 1. The black squares and circles show 
At corresponds to the case of the transfer fidelity of an excita- 
tion and the average fidelity with the system subjects to dissi- 
pative environment, while the red crosses and plus signs show 
At corresponds to those with the system subjects to dephasing 
environment. Moreover, the average fidelity was plotted with 
the phase shift being corrected, and the inset is a schematic 
picture (with the system subjects to dissipative environment) 
annotating the meaning of the FWHM. 



a quantum state with the highest fidelity 2/3 |46j. The 
typical plots of 7c versus N with scaling constant A = 1 
is shown in Fig. 3, where we have assumed that the phase 
shift caused by the factor cos a is eliminated by choosing 
cos a = 1 for the cases of iV = 4r — 1 and N = 2r. From 
this figure one can see that 7c decreases with the increase 
of N, this puts another constraint on this spin channel for 
long distance quantum communication because for large 
A'', one needs to reduce the mutual interaction between 
the system of interest and its surrounding environment in 
order to obtain a reliable transfer fidelity, and in general, 
this is not an easy task in real experiments. 

From the aforementioned analysis one can see that in- 
creasing the decoherence rate 7 or the spin chain length N 



always introduces stringent constraints on the spin chan- 
nel for long distance quantum communication. In the fol- 
lowing, we discuss another serious constraint imposed by 
the dissipative and dephasing environments on the spin 
channel. For this purpose, we examine full width at half 
maximum (FWHM) of the first oscillation describing the 
dynamics of the fidelity of an excitation as well as the 
average fidelity and we denote it as At = \t2 — ti\, with 
fih) = f{t2) = frnax/2 and F(ti) = ^(is) = Fmax/2. For 
the case of dissipative environment, from Eq. (26) one can 
obtain the following two relations 

e-7tsin2(w-i)At = /„,„,/2, 



e~''*/2 sin^^i At \/3F.,nax - 3 + cos2 a - cos a, (27) 

where for transfer fidelity of an excitation, (^1,^2) can be 
obtained by solving the nonlinear equation expressed in 
the first line of Eq. (27) , and for the average fidelity, {ti , t2) 
can be obtained by solving the second equation expressed 
in Eq. (27). The maximum value of fmax and F^ax can 
be obtained directly from Eq. (26). 

For the case of dephasing environment, however, one 
cannot obtain a result similar to that of Eq. (27), thus we 
resort to numerical methods. The exemplified plots are 
shown in Fig. 4, where for the case of average fidelity, we 
have chosen cos a = 1. From this figure one can see that 
the curves exhibit nearly the same behaviors for both of 
the decoherence scenarios. With the increase of the spin 
channel length N, the FWHM At decreases one by one, 
which implies that the peak of the wave becomes steeper 
and narrow with the increase of N. This puts a stringent 
limitation on the time when the measurement must be 
performed by the receiver Bob, because for large N, when 
the transmitted state arrives at the destination qubits. 
Bob should make measurements on the qubits at the in- 
stant time of tc as quickly as possible, otherwise, even a 
minor deviation of measurement time from tc would lead 
to a great inaccuracy to the transmitted state. 

We now turn to investigate quantum state transfer in 
another spin channel proposed by Shi et al. [6] , the Hamil- 
tonian of the system has the similar form as that expressed 
in Eq. (19), with however, the coupling strength between 
neighboring spins given by 



A-\/ ■n{N — n) [n £ even), 

Xy/ (n + 2k){N - n + 2k) {n G odd). 



(28) 



Here k = {0,1,2,...}, and A is still a scaling constant. 
This Hamiltonian is the same as that of Eq. (19) when 
k = 0. For other cases of k and N e even, it can be 
used to achieve perfect state transfer from the first node 
to the last node at time to — tt/2X and periodically re- 
turns there at regular time intervals of tt/A in the ab- 
sence of decoherence environment (for the special case of 
N — 2, perfect state transfer can also be achieved at time 
t = TT/{4:k + 2)X + TTr/{2k + l)X, r = {0, 1, 2, . . .}). Here we 
discuss how the dissipative and dephasing environments 
affect quantum state transfer in this spin channel. Since 
the case of fc = has been discussed in the above section. 
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here we concentrated only on dynamics of state transfer 
fidelity at time to — {2k — l)7r/2A for the case of fc ^ 1. 
The eigenvalues and the eigenvectors of the correspond- 
ing Hamiltonian have been obtained in Ref. [5], and this 
enables an analytical analysis for the case of dissipative 
environment, while for dephasing environment, we also 
resort to numerical analysis. Our results revealed that for 
both of the decoherence modes, the transfer fidelity of an 
excitation and the average fidelity at time to for fc ^ 1 
show completely the same behaviors as that of k — and 
TV € even, thus the region in which this spin channel is su- 
perior to any classical communication channel is also the 
same as that of A; = and N G even. However, for any 
fixed N e even, the FWHM At becomes narrower and 
narrower with the increase of k, thus by the same argu- 
ment as made in the previous section, we suggest that in 
practice, one can choose A: = in order to minimize the 
deviations from the transmitted state introduced by the 
measurement process. 

Moreover, for the special case of = 2 at arbitrary 
evolution time t, the transfer fidelity of an excitation and 
the average fidelity with the system subjects to dissipa- 
tive environment can be obtained analytically as f{t) = 
e-T*sin2[(2fc + l)Ai] and ^ e'^Hin^[{2k + l)\t]/6 + 
1/2. For the case of the system subjects to dephasing envi- 
ronment, again due to the fact that yl(2)] 7^ 0, one 
cannot obtain an analytical result of the state transfer 
fidelity. 



4 Entanglement distribution and creation in 
decoherence spin channels 

From the above discussion one can see that both the dis- 
sipative and dephasing environments have severely detri- 
mental effects on fidelity of quantum state transfer. How 
these decoherence modes affect other dynamical processes 
of a quantum system. As an answer to this question, in this 
section we investigate decoherence effects on entanglement 
distribution and creation. We first see entanglement distri- 
bution between two distant parties through the dissipative 
and dephasing spin channels. To assess the extent to which 
these decoherence environments affect entanglement dis- 
tribution, we assume the entangled state \ip) = (|01) -I- 
|10))/-\/2 is initially prepared between a non- interacting 
qubit NI and the first qubit A of the chain, then after some 
time t, the entanglement may be established between NI 
and the A'^th qubit B. The overall Hamiltonian of the sys- 
tem can be written as Hover = I2 ® H {I2 denotes the 
2x2 identity matrix). To assess the amount of the pair- 
wise entanglement at different instants of time, we adopt 
the concurrence, a function introduced by Wootters [47 , 
equals to 1 for maximally entangled states and zero for 
separable states, defined as C = max{0, Ai — A2 — A3 — A4}, 
where Ai {i — 1, 2, 3, 4) are the square roots of the eigen- 
values of the matrix R = Pmnio'm ^n)Pmn(^m ® ^n) 

{pmn is the two-qubit reduced density matrix). Then af- 
ter a similar numerical analysis as performed in Section 3, 
we obtained entanglement dynamics of the system. The re- 



sults revealed two remarkable features (see Fig. 5). First, 
for both of the decoherence modes, the pairwise entan- 
glement between the non-interacting qubit NI and the 
last qubit B measured by the concurrence Cni.bI^) dis- 
play completely the same behaviors, i.e., they both be- 
have as suppressed fluctuations as the time t evolves, and 
the entanglement may experience sudden death 148) dur- 
ing certain intervals of the evolution time (see the inset in 
the top panel of Fig. 5). Second, for any fixed decoherence 
rate 7, the maximum value of Cni.b(0 decreases smoothly 
with the increase of the spin chain length N, and the crit- 
ical time tc at which the concurrence Cm.B{t) attains its 
maximum value still does not locate at to = 7r/2A (at 
this point, the concurrence can be obtained analytically 
as CNi,B(io) = e~^*'>, which is independent of N). tc in- 
creases as N increases, and when iV — >■ 00, approaches 
its asymptotic value 7r/2A, which is the same as that of 
to- This implies that even for very large N, one can still 
obtain substantial concurrence Cm,Bitc) — e~^'^l'^'^ (it is 
not difficult to understand this phenomenon, because the 
scaling for large iV also involves increasing the coupling 
strengths between the neighboring spins for we have cho- 
sen Jn,n+\ = \^n{N - n)). 

Next we turn to investigate how the dissipative and de- 
phasing environments affect the creation of entanglement 
in spin networks. For this purpose, we consider the multi- 
arm structure AI{Ni, N2, Na) of the XX spin chain Hamil- 
tonian (19) with the addition of the exchange couplings 
between the hub site and its nearest-neighbor output sites 
satisfy the branching rule J,0,49 . Here iVi and A'2 denote 
the number of sites in the input and output arms, respec- 
tively, and Na denotes the number of output arms (see the 
inset in the top panel of Fig. 6). It has been shown that 
in the absence of decoherence environment, this structure 
can be employed to create multi-qubit entangled W state 
[lOllSO] at the ends of the outgoing arms. When the system 
is subjected to the decoherence environment, however, the 
case may be different. 

When the quantum system subjects to dissipative en- 
vironment, our numerical results revealed that the concur- 
rence C (here C measures amount of the created pairwise 
entanglement between arbitrary two end nodes of the out- 
put arms) is determined by two numbers iV = TVi -I- iV2 -I- 1 
and Na- For any fixed Na and decoherence rate 7, the 
maximum value of the concurrence decreases with the in- 
crease of N- This is different from that of the ideal case 
(i.e., 7 = 0), for which the concurrence has no relation 
with TV [TUifiU] . Moreover, similar to the case of entangle- 
ment distribution in a spin chain, the critical time tc at 
which the concurrence attains its maximum value is also 
slightly different from to = 7r/2A for small values of N (see 
the top panel of Fig. 6), while tc is nearly independent of 
Na- In the limit of TV 00, however, tc ^ to — 7r/2A, 
at which the concurrence can be obtained explicitly as 
C{tc) - 2e"('^'^/2^V^A. 

When the system subjects to dephasing environment, 
however, our numerical results revealed that the concur- 
rence C{t) is dependent on iVi, N2 and Na, which is dif- 
ferent from that of the dissipative environment. For fixed 
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Fig. 5. (Color online) Rescaled critical time tc/ir and the 
maximum concurrence CNi,B(tc) versus with the system 
subjects to dissipative or dephasing environment. For every 
plot, the diamonds from top to bottom correspond to 7 = 
0.1,0.2,0.3,0.4,0.5. The inset in the top panel shows dynam- 
ics of the concurrence Cni.b versus t/vr. The other parameter 
values are given by 7 = 0.1 and A = 1. 



N — Ni + N2 + 1, the maximum value of C{t) decreases 
with the increase of N2, thus in order to minimize the 
detrimental effects of dephasing environment on entangle- 
ment creation, one can design the spin structure with the 
number of sites in each output arm as N2 — 1. Moreover, 
as can be seen from Fig. 6, when the system subjects to 
dephasing environment, the critical time at which the 
concurrence attains its maximum approaches the asymp- 
totic value 7r/2A more quickly than that for the dissipative 
environment, and the detrimental effects of the dissipative 
environment seems to be more severe than that of the 
dephasing environment under the condition of the same 
decoherence rate 7. 



5 Summary and discussion 

In the present paper, by exact diagonalization of a gener- 
alized superoperator for solving the master equation, we 
investigated quantum state transfer, as well as entangle- 
ment distribution and creation through the dissipative and 
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Fig. 6. (Color online) Rescaled critical time tc/n and the max- 
imum concurrence C{tc) versus A^i with the system subjects 
to dissipative (black diamonds) and dephasing (red squares) 
environments. The other parameters for the plots are N2 = 1, 
Na = 3, j — 0.3 and A = 1. The inset in the top panel is an 
exemplified multiarm structure of the spin network with the 
solid red circle denotes the hub, and the number of sites in the 
input and output arms are given by A^'i —4, N2 = 2, while the 
number of output arms is Na = 3. 



dephasing spin channels, aimed at revealing the extent to 
which these decoherence environments affect various QIP 
tasks. We focused on several interaction-modulated spin 
networks which may serve as perfect spin channels in the 
absence of decoherence, and gave an analytical analysis for 
the case of dissipative environment, and numerical anal- 
ysis for the case of dephasing environment. The results 
revealed three general conclusions. The first is that these 
two decoherence environments always lead to suppression 
of the state transfer fidelity, as well as the amount of the 
distributed and created pair wise entanglement, particu- 
larly, this suppression increases with the increase of the 
decoherence rate 7. The second is that the detrimental ef- 
fects of the dissipative environment are more severe than 
that of the dephasing environment under the condition of 
the same decoherence rate 7. The third is that the criti- 
cal time tc at which the transfer fidelity and the concur- 
rence attain their maxima increases with the increase of 
N, and when A'^ — cxi, it approaches the asymptotic value 
to = 7r/2A. 

For the case of state transfer, the results also revealed 
that for fixed 7, the transfer fidelity of an excitation /(to) 
decreases with the increase of while the average fidelity 
F{to) displays three different patterns corresponding to 
the chain length N ^ Ar + 1, N = Ar - I and N = 2r, 
and only when A'^ = 4r -f- 1 can the spin channel superior 
to all of its classical counterpart if there are no additional 
operations being performed. Besides this, when the system 
subjects to dissipative environment, /(<o) is independent 
of A^, while for each pattern A^, F{to) is independent of 
r. When the system subjects to dephasing environment, 
however, both f{to) and F{to) decrease with the increase 
of N and r. Moreover, the critical decoherence rate 7c 
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before which the decoherence spin channel is superior to 
its classical counterpart is independent of N for dissipative 
environment, and it decreases with the increase of N for 
dephasing environment. 

For entanglement distribution and creation, our results 
revealed that the distributed entanglement shows com- 
pletely the same dynamical behaviors when the system 
subjects to dissipative or dephasing environment, and for 
fixed 7, the maximum value of the concurrence Cni,b de- 
creases as N increases, and when iV — >■ oo, CNi,B(ic) = 
g-77r/2A^ The created entanglement also decreases as N 
increases, and for the case of dissipative environment, the 
maximum concurrence C approaches its asymptotic value 
C{tc) = 2e-('y'^/2^V^A in the limit of N ^ oo, while for 
the case of dephasing environment, its asymptotic value 
can only be obtained numerically. 

From the above analysis it can be seen that both the 
dissipative and dephasing environments have severe detri- 
mental effects on transport capacity of a spin channel, thus 
a question naturally arises at this stage is how to stabilize 
a quantum system against this unwanted phenomenon, 
or in other words, how to minimize, delay or even elim- 
inate the detrimental effects of decoherence. In a recent 
work [51) . we have shown that for some specific environ- 
ments, it is possible to obtain a near perfect state transfer 
by adjusting the coupling strengths between neighboring 
qubits of a spin network, however, for decoherence modes 
considered here, this procedure cannot be applied, thus 
one must resort to other techniques in order to get a high 
efficient communication. Physically, one can increase the 
transfer fidelity to some degree by increasing the coupling 
strength between neighboring qubits of the spin network, 
or equivalently, by increasing the scaling constant A (this 
statement can be proved true by the explicit form of the 
operator A, from which one can obtain that increasing A 
will decrease the period of the suppressed oscillation and 
thus enlarges the maximum value of the state transfer fi- 
delity). Particularly, in the limit of A — > (X), the transfer 
fidelity will approaches unity. However, the realization of A 
large enough to obtain high efficient communication might 
be a difficult experimental task. 

Finally, we would like to emphasize that although there 
are works showing that efficiency of the quantum state 
transfer can sometimes be enhanced to some extent by 
decoherence [521 I531 I54I, finding ways to minimize or even 
eliminate the detrimental effects of decoherence is still the 
prerequisites and challenging task in the practical realiza- 
tion of quantum computer. 
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